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Abstract: Given any compact hyperbolic surface M, and a closed 
geodesic on M, we construct of a sequence of quasimodes on M whose 
microlocal lifts concentrate positive mass on the geodesic. Thus, the 
Quantum Unique Ergodicity (QUE) property does not hold for these 
quasimodes. This is analogous to a construction of Faure-Nonnenmacher- 
De Bievre in the context of quantized cat maps, and lends credence to 
the suggestion that large multiplicities play a role in the known failure 
of QUE for certain "toy models" of quantum chaos. We moreover con- 
jecture a precise threshold for the order of quasimodes needed for QUE 
to hold — the result of the present paper shows that this conjecture, if 
true, is sharp. 



1. Introduction 

An important problem of "quantum chaos" is to understand the 
relationship between a dynamical system and its high-energy quan- 
tum mechanical models. The "correspondence principle" dictates that 
quantum mechanics should replicate classical mechanics in the semi- 
classical limit ft — > 0, and this is generally the case when the classical 
system is completely integrable. In situations where the dynamics are 
more complicated, however, the standard analytic techniques are only 
valid for "short" times, up to the scale of the Ehrenfest time C\ log ft |, 
where C is a constant depending on the system. 

Consider the geodesic flow on a compact surface M = V\M of con- 
stant negative curvature. It is well known that this dynamical system 
is "chaotic" — eg., mixing (with respect to Liouville measure on S*M), 
Anosov, etc. On the other hand, the high-energy spectral data for 
such surfaces is extremely mysterious (see eg. [Sar03j ). The Quantum 
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Unique Ergodicity (QUE) Conjecture states that high energy eigen- 
functions become equidistributed as the eigenvalue tends to oo — cor- 
responding to the semiclassical limit h — > 0; precisely, that the distri- 
butions 

f M <t) :feC 00 (S*M)^(Op(f) ( j ) , ( p) 

converge in the weak-* topology to the Liouville measure on S*M as 
the eigenvalue of tends to oo. Any weak-* limit point of the /i^ is 
called a quantum limit . 

It is known by |Sni74 IZel87t ICdV85] that any quantum limit is a 



(positive) measure invariant under the geodesic flow; if we normalize 
||0|| 2 = 1, then they are probability measures. They further show that 
almost all /i^ become equidistributed, in the following sense: for any 
choice of orthonormal basis of L 2 (M) consisting of Laplace eigenfunc- 
tions, there exists a zero-density exceptional set of basis eigenfunctions 
such that the remaining eigenfunctions satisfy QUE. This property is 
known as Quantum Ergodicity, and holds in great generality — it 
depends only on the ergodicity of the flow. QUE then asks that there 
be no exceptional subsequences. 

In contrast, there are "toy models" of quantum chaos that demon- 
strate non-QUE behavior. For example, hyperbolic linear maps A G 
SX(2,Z) of T 2 — also known as "cat maps" — can be quantized, and 
shown to satisfy Quantum Ergodicity [BDB96J . Yet not every quantum 
limit need be Lebesgue; examples are constructed by Faure-Nonnenmacher- 
De Bievre [FNDB03J of quantum limits that are half-Lebesgue and half 
atomic (i.e., half of its mass is supported on a finite periodic orbit). 

It is suspected [Sarllj that large spectral multiplicities play a role in 
this phenomenon. For the cat maps, the quantum propagator has finite 
order, given by ordjv^) where A G SX(2,Z) is the matrix generating 
the classical dynamics, N = (2irh)^ 1 G Z is the inverse Planck's con- 
stant, and ordAr(v4) is the order of the matrix A modulo N. It is known 
that ord7v(v4) > 21og A N + 0(1), where A > 1 is the absolute value of 
the large eigenvalue of A. The examples of |FNDB03] rely crucially 
on minimal periods of the propagator; the sequence of eigenvectors is 
taken from values of N for which ordAr(A) saturates this lower bound. 
Since the propagator is unitary, the eigenvalues must be ordAr(A)-roots 
of unity; the Hilbert space of states from which these eigenvectors are 
taken is iV-dimensional, and so such short periods imply massive de- 
generacies in the spectrum of size about 

N 1/2 
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In fact, Bourgain shows |Bou07] (improving on earlier results of Kurlberg- 
Rudnick |KR01j ) that QUE holds whenever ord N (A) > t > h~ e . 
However, there is still a significant gap between the logarithmic scale 
0(1/| log ft |) of the non-QUE examples and Bourgain's QUE bounds. 

In fact, we conjecture the correct threshold to be just beyond the 
logarithmic scale of the known counterexamples: 

Conjecture 1. Let n be a quantum limit for the cat map A, such that 
the orders of the propagators satisfy 

hminf 0r ^ A K C 

7V->oo log N 

Then the entropy of fi is bounded below by 

1 



In particular, if 



hip) > h max - 



,. . r ord N (A) 
hmmt — = oo 

JV-voo log N 



then fi is Lebesgue measure. 

We now return to the geodesic flow on a compact hyperbolic surface 
M. We parametrize the spectrum of the Laplacian on M by A = \ + r 2 , 
and here the frequency r — > oo plays the role of the inverse Planck's 
constant in the semiclassical limit. It is conjectured |IS95] that the 
multiplicities are bounded by O e (r e ), but this problem is well out of 
the reach of current technology; in any case we do not have control 
over different multiplicity scales here as in the case of cat maps. 

However, one can artificially introduce "degeneracies" by consid- 
ering quasimodes, or approximate eigenfunctions. Define an u(r)- 
quasimode with approximate parameter r to be a function ip 
satisfying 

||(A + (i + r 2 ))^|| 2 <ra;(r)||^|| 2 

The factor of r in our definition comes from the fact that r is essentially 
the square-root of the Laplace eigenvalue. Intuitively, we think of ip as 
being localized spectrally near the window [r — u(r), r + oo(r)]. 

The main term in Weyl's Law (see eg. [Ber77]) says that the as- 
ymptotic density of eigenfunctions near spectral parameter r is pro- 
portional to r (with a constant depending on the area of the surface) . 
Though controlling the error term is a very difficult problem — indeed, 
this is precisely the problem of bounding multiplicities — it is known 
|Ber77] that this approximation is valid for "large logarithmic windows" 
cj(r) > Kj logr, where again K depends on the surface, and is believed 
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to be valid for far smaller windows; eg., u(r) = r~ e . In any case, we 
shall be interested in this paper with windows of size u(r) > 1/logr, 
where we certainly expect the error in Weyl's Law to be negligible. By 
varying the size u(r) of the windows, we can analyze its role in QUE 
phenomena. 

Conjecture 2. Let M be a compact hyperbolic surface, and {ipj} a 
sequence of o (■^-^) -quasimodes. Then {ipj} satisfies QUE. 



This conjecture is a slight strengthening of the QUE Conjecture 
[RS94] for the case of Riemann surfaces, which essentially claims (in 
a quantitative way) that the small spectral multiplicities are respon- 
sible for QUE in this setting. It is also important to observe that 
the windows u(r) = o(l/logr) are just beyond what can be analyzed 
at present: as remarked by Sarnak [Sarllj . any proof of QUE is ex- 
pected to address the multiplicity issue (though perhaps indirectly). 
Conjecture |2] suggests that QUE is on par with a multiplicity bound 
of o(l/ log r), whereas the current best known bound due to Berard 
|Ber77] is 0(l/logr). 

In this paper, we study "logarithmic-scale" quasimodes, by which we 

mean ( j-^— ) -quasimodes for some constant e > 0. Precisely, we prove 



Theorem 1. Let M be a hyperbolic Riemann surface, and £ C S*M a 
closed geodesic. Then for any e > 0, there exists S > and a sequence 

of ^1^7 j -quasimodes on M whose microlocal lifts do not equidistribute; 

indeed, they concentrate mass > 5 on the geodesic £. 

Thus, if true, Conjecture [2] is sharp. This is reminiscent of — and, 
indeed, the argument inspired by — the constructions of |FNDB03] for 
cat maps. 

We also remark that the case of arithmetic joint quasimodes is radi- 
cally different, as it was shown in |BL12] that any sequence of joint 
o(l)-quasimodes satisfies QUE, even though this includes spaces of 
much larger dimension o(r) (with arbitrarily slow decay). This is due 
to the additional rigid structure imposed by the arithmetic symmetries, 
which are not respected by the quasimodes constructed here in Theo- 
rem [T] (they a fortiori cannot be o(l)-quasimodes of a Hecke operator). 

The main part of the argument is constructing Cj log r-quasimodes 
that work, for some large C; moving from a large coefficient to a small 
e is achieved by means of a (somewhat cheap) argument to dilute out 
the Cj logr-quasimode, given in section |5j The construction of these 
quasimodes for a sufficiently large C is the subject of section [2j and is 
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based on a construction from [Brol2j . In section [3] we compute the sta- 
tionary phase asymptotic which is the analytic heart of the argument, 
and then apply it in section H] to prove that our examples are indeed 

O ^1^7 j -quasimodes that localize a positive proportion of their mass 

on the geodesic £. 

Though the technical approach is slightly altered for simplicity, in- 
tuitively the main idea is to begin with a radial 1-quasimode localized 
near (more precisely, localized near radial vectors in a neighbor- 

hood of p), and average over propagation times to improve the quality 
of the quasimode to order C / logr. To get localization near £, we then 
average these spherical quasimodes over a piece of the stable horocycle 
through p, which enhances the quasimodes along £ through construc- 
tive interference, while canceling out the mass away from £ through 
destructive interference. Naturally, then, the main analytic tool for ex- 
ploiting the interferences will be the method of stationary phase, which 
will yield asymptotics for the relevant integrals defining the microlocal 
lifts of the quasimodes. 

Acknowledgements. The question of proving the existence of such 
quasimodes was first posed to the author by Peter Sarnak. We also 
thank Lior Silberman and Elon Lindenstrauss for helpful discussions. 



We begin with a generalization of the discussion in |Brol2j . where 
we exhibited spherical kernels whose microlocal lifts localize near radial 
vectors. Intuitively, we wish to take a localized kernel and improve the 
order of the quasimodes by averaging over propagation times with an 
operator like 



which cancels out components of spectral parameter away from r. 
This is reminiscent of — and largely inspired by — the construction 
in [FNDB03J of scarred eigenstates for quantized cat maps. To get an 

O ( t — ) -quasimode, we would have to average up to time T > logr. 



From a technical perspective, though, it is easier to proceed in the 
following manner — which is essentially just a smoothed version of this 
averaging procedure. Recall the Selberg/Harish-Chandra transform 
for a spherical kernel A; on EI (see eg. |Iwa02l Chapter 1.8]), giving the 
eigenvalue h(s)<f> 8 = k * (fi s for a Laplace eigenfunction <f> of spectral 



2. Construction of the Quasimodes 
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parameter s: 

h(s) = I j su g{t)dt 

t 
2 

1 f°° dQ{u 



g(t) = 2Q ( sinh 

;i) k{ u ) = 



7T J u ^UJ-U 

The coordinate u(z, w) = smh. 2 (d(z, w)/2) is often more convenient for 
calculations than the actual distance (du is the radial volume measure 
on H). What is most important for our purposes is that whenever g 
is compactly supported in the interval [— T, T], the kernel k will be 
supported in the ball of radius T in H. 

We can write such a kernel k as a (left-i^- invariant) function on H, 
and use Helgason's Fourier inversion [Hcl81J to write 

fe(^)= / / e {is+ ^ <z ' b> k(s,b)stanh(7rs)dsdb 

Js=0 JB 

where b G B runs over the boundary S 1 of the disc model for H, 
and < z,b > represents the (signed) distance to the origin o from the 
horocycle through the point z G EI tangent to b G B. Since each 
plane wave e < -~* s+ 2) < ' b> is an eigenfunction of spectral parameter s, 
the Fourier transform 

k(s,b) = [ e { - is+ ^ b> k{z)dz 
Jw 

= h(s)e { - ts+L z )<0 > b> 
= h(s) 

so that 



k(z) 



J (^j e {is+1 i )<z > b> db^j /i(s)stanh(vrs)c/s 



It will be more convenient to write as in |Zel87] 

e{is+ i )<z , b>db = e (is-±)<z,b> d9 = e {is-\) V {g.k e ) de 

where tp{g) is again the signed distance from the origin to the horocycle 
through g G PSL(2, R), and k e parametrizes the SO (2) fibre gK. Since 
if is left i^-invariant and right N- invariant, it is sometimes convenient 
to use KAN coordinates to write 

9 ~[sm9 cos9 J I e~*/ 2 J\0 1 J [c d 
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so that this distance is given by <p(g) = t = log (a 2 + c 2 ). 

Fix an orthonormal basis {<fti} of L 2 (M) consisting of Laplace eigen- 
functions, which we can take to be real-valued for simplicity. Each 
eigenfunction generates, under right translations, an irreducible rep- 
resentation Vi = {^(xg- 1 ) : g e PSL(2,R)} of PSL(2,R), which to- 
gether span a dense subspace of L 2 (S*M). 

We distinguish the pairwise orthogonal weight spaces A 2n in each 
representation, consisting of those functions satisfying f(gkg) = e 2m0 f(g) 

for all kg = ( C ° S f ) E K and g G S*M. The weight spaces to- 



sin 6 cos 9 

gether span a dense subspace of V/. Each weight space is one-dimensional 
in Vi, spanned by 02« where 

4° = <j>ieA 



(i n + ±+n)0£+2 = E+cli^ 



[in + - - n)^l_ 2 = E-^ n 



1 

2 
1 

2 

Here E + and are the raising and lowering operators, first- 

(\ i \ 

order differential operators corresponding to I ^ ^ 1 e sl(2,M) and 
1 -A 

^ 1 ) ^ ^(2, R) in the Lie algebra. The normalized pseudodiffer- 
ential operators 

• 02n ^ 02n-2 

are unitary and left-invariant, and each 02n is a unit vector. We define 
the distribution 



*2= E 



4 (0 

; 2n 



and extend this definition by linearity to = Yl™=-oc = S^L-oo 
for quasimodes ip, where each ip 2n = R n i/i- 

2.1. Construction of the Microlocal Lifts. The following can be 
found in |Brol2j . based on the arguments of [LinOl] . We set 

W) = (Op(f)M) ■= (/*oo,V>> = ( / V ^2n,^0 

\ re=-JV 

according to the pseudo-differential calculus of [Zel87], which clearly re- 
stricts to the measure \ip(z)\ 2 dz when applied to i-T-invariant functions 
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/ G C°°(M), by orthogonality of the weight spaces. Note moreover 
that this limit is purely formal for i^-fmite /, and since these K- finite 
functions are dense in the space of smooth functions, we can restrict 
our attention to these. 

Lemma 1. Let i/j be a linear combination of eigenf unctions with spec- 
tral parameter in [r — l,r + 1], with \ \1pW2 = 1, and set 

I 3L v-^ L — \n\ . 

\n\<L 

Then for any smooth f G C°° (S* M) , we have 

l4f) = (fy i y) + O f (L~ 1 + Lr- 1 ) 

In particular, \I$(f) — \&) | —> if we set L = \r a \ for some 
< a < 1. 

Proof: See [LmOTl IBLT21 IBrol2j . 

Note that the prefactor 2 l^+i ~ ^/f r ~^ 4 ^ s shnply an L 2 -normalization 
of the Fejer coefficients . 

Corollary 1. Let {ipj} be a sequence of o(l)-quasimodes with approx- 
imate spectral parameters rj, normalized in L 2 (M), and L = rj a for 
some fixed < a < 1 as above. Then for all f G C°°(S*M) we have 

IW/)-(/*i,*i>|->0 

Proof: As in the proof of Lemma [1] from |Lin01l IBrol2j , we may 

take / G ^|n|<7V ^ 2n ^° ^ e ^-&mte, since these span a dense subspace 
of C°°(S*M). The Lemma shows that the statement holds for the 
projection of ipj to the space spanned by eigenfunctions of spectral 
parameter in [rj — 1, Tj + 1]; it is sufficient to check that the contribution 
of other spectral components is negligible. We decompose ip = ip- in + 
■0out, an d the corresponding \1/ = + \I/ out , where the "in" component 
consists of all spectral components inside the interval [r — l,r + 1], 
and the "out" component is the orthogonal complement consisting of 
spectral components outside this range. 
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We have 

-(/* in ,* ill ) < 



< 
< 



(/*,*aut) 



(Z^out, 



/ 



I* 



out 2 



2 ' H^out || 2 

HV'outlb 



out 1 1 2 



I* 



in 1 1 2 > 



\n\<N 



out 



/ (^out) 2 „ , "0u 

|n|<AT 



< (2iV + 1) 

HVwIh 



(IMMWUIa + HV'. 



out 2 ' 



% 2 



recalling that each | \ip2n\ 1 2 — 1 1^| I2 by unitarity of the raising and lower- 
ing operators R and -R -1 , and so it remains to show that | iVwIh = o(l). 
But note that for any spectral component <fii with spectral parameter 
outside [r — 1, r + 1], the defect ||(A — (7 +r 2 ))0/| | 2 > r| |0/|| 2 , and thus 

ll(A-(^ + r 2 ))^l|2>||^out||2-r 

so that the o(l)-quasimode condition necessitates H^outlh = o(l), and 
we get the desired asymptotic. □ 

Thus, for any given sequence {ipj} of o(l)-quasimodes (and in par- 
ticular, for sequences of 0(1/ logr)-quasimodes), we have constructed 
a sequence such that the microlocal lifts \^/j\ 2 dVol are asymp- 

totically equivalent to the distributions 1^.. It is these measures that 
we wish to study. 



2.2. The kernels k^j. Propagating spherical kernels over long times 
will equidistribute on S*M, but by averaging over a stable neighbor- 
hood of the geodesic £, we will obtain quasimodes that localize a posi- 
tive proportion of their mass on £. The intuitive picture is that averag- 
ing along a stable horocycle causes constructive interference at vectors 
pointing along £, enhancing the mass near the geodesic, while caus- 
ing destructive interference away from £. Naturally, the basis for our 
analysis in the coming sections will be the method of stationary phase. 

Thus, we pick a point p G £ on our geodesic, and choose a covering 
EI — > M such that the origin % e HI is mapped to p, and the imaginary 
= is mapped onto the geodesic £. Let {rj} C • 2ttZ be a 
sequence of resonant frequencies for the length /(£) of the geodesic £; 
this condition on the approximate spectral parameters will ensure that 
our quasimodes will not self-interfere as they wrap around £, and will 
be used in the calculations of Proposition [3j 
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We then take the spherical kernel given, via §Bj, by the Fourier pair 

r, , COsh zfc cosh it 3 
his) = — - — — p- 

COSh -rr- + COSh 

m = k CM(tr > ) 



2 J cosh (ttKjO 

for fCj := 2k§7"- We then choose a smooth, positive, even cutoff func- 



tion x supported in [—1,1], which we may assume is identically 1 on 
[-1 II 

L 2' 2-1' 



1 1 1 and replace 



/i( s ) = h(s)*x 



log r j 



C 

so that g is supported in the interval of radius logr^. We denote 
the resulting kernel by kj, and its micro local lift — as constructed in 
Lemma [T] — by Kj] they are supported in the ball of radius C~ x log rj 
in H and S*M respectively. Note as in |Brol2] that the distribution 

oo „ 

R n kj(9)= e {ls -^ M9 - ke) deh(s)stanh(TTs)ds 

n=-oo s 

by |Zel87] , with kg = ^ cos0 y ' an< ^ SO we ma ^ wr ^ e K i as a 
(right-)convolution of Yl^=-oo with the L 2 -normalized Fejer kernel 
of order L on the 5*0(2)-fibre K: 

[ [ e (is ~^ M9 - k(,) F L (29)deh(s)staBh(TTs)ds 

Js J 8=0 



where F L (26) := E W <x «» (2n0). We will take L 

[rj ^] as a convenientEl choice for L. 



We note for later use the estimate 



< r" 1 



(2) / S - 1 /i(s) tanhvrscfe < — - - ^_ , , 

Js °Cj O 

and remark, as in [Brol2j . that since x £ C°° is smooth, combined with 
the rapid decay of /i away from Tj, the function /i(s) is also rapidly 
decaying away from Tj\ the estimate h(s) < \rj — s\~ 3 will suffice here. 



1 It should be emphasized that here and throughout we make no attempt to 

obtain the sharpest possible bounds, and thus allow ourselves to be wasteful with 

1 IC 

powers of r- ' , at the expense of taking a larger constant C than necessary. 
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We now define 

/ ' n 



hA z ) = J k \)' z ) x ^ dn 

to be the convolution of kj with the smooth cutoff \i along the stable 
subgroup N. Since this convolution acts in kj on the left, it commutes 
with the left-invariant operators R and defining our microlocal lift, 
and so K^j is given by 



K tM = I x(n)K(n.g)dn 

x(n)F L (2a)e {is -^ {n - 9 - ka) dadn) h(s)s tanh(vrs)rfs 




n J aGT 



We also define 



z 



7er 



to be the respective projections to M and S*M. Note that since k^j 
and K^j have compact support, the sums are finite; and that as usual 
||%ilU 2 (M) — ||^,il|i/ 2 (s*Af)- Naturally these projections do not com- 
mute with the averaging over the stable horocycle that defines k^j and 
k^j, but it will still be easier for us to approximate pointwise, and 
then estimate the sum over V to project back to S*M; in reality, the 
vast majority of the values of n^j in the sum over T will be negligible. 
We will show that: 

• For any neighborhood U of the geodesic £ in S*M, we have 
j v \R u \ 2 dVol > l/log7> 



• The full L 2 -norm satisfies < 1/log 

• The function k^j is an O ^j^^-j-quasimode on M. 

Here and throughout, we allow all implied constants to depend on C . 
Together these prove Theorem [I] for the case of a large constant; in 

the final section we show how to bootstrap the argument to (j^r. 

quasimodes, for arbitrarily small constants e > 0. 



3. A Stationary Phase Estimate 
We now perform the main stationary phase approximation that is 



the crux of our analysis. We write ni s ] for the s-spectral component of 
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K^j] that is to say 



J s 

K$(x,t,0) = [ [ x(x-n)F L (2\9-a\)e {ts -^ {nAa) dadn 



n J a£T 



where we denote 



ip{n,t,a) = pll 1 n \( 6 e _ t/2 



cos a — sin a 
sin a cos a 



' e t//2 cos a + ne i//2 sina * 

V I P ~t/2 ■ 



e ' sin a * 

2„ i • /o„A i „-t/„2 i i\ • 2 



(3) = log(e cos a + nsin(2a) + e (n + 1) sin 2 a) 



in NAK coordinates. We will perform the stationary phase analysis 

(s) 

in the integrals over n and a, estimating j pointwise, and then deal 
with the integral over the spectrum afterwards. 

Proposition 1. Let (x,t,9) e S*M with \t\ < C^hgrj, and set L = 
[r^°^ c \ . Then 



K W( X) t, 9) = 7rs- 1 x(x)e t/2 [e lst F L (26) + e- lst F L (ix - 29)] + O(r] 00/C ' s~ 2 ) 

Remark: Along the geodesic, the e t//2 term represents the growth due 
to constructive interference in the stable direction. This exponential 
term will play a central role in the mass concentration of Proposition [3j 

Proof: The proof is an application of the method of stationary 
phase (see eg. |Hor03t IZwol2j ). Note that is defined pointwise 
via integrals over N and K; the integral over N takes place on the 
left, while the integral over K takes place on the right, so that these 
translations commute. The sole critical points of the phase func- 
tion ip are at the origin (n,a) = (0,0) and at (n, a) = (0, 7r/2); to 
see this it is actually easiest to work in .fTAiV-coordinates, to take 
advantage of the fact that (p is left-i^ and right- N invariant. Let 

{cos 9 — sm9\ (e l l 2 \ (\ u\ nnr/n m\ • i 

3 = [sin9 cos* ) { e-/ 2 ){ l) e ^L(2,R); we wish to 

calculate the critical points of ip for translation along N on the left, 
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and K on the right. Thus we compute 



d_ 

dn 



(p(n.g) 



n=0 



d_ 

dn 
d_ 

dn 
d_ 

dn 

d 



1 n 
1 

1 n 
1 



cos 9 — sin 9 

sin 9 cos 9 

cos 9 — sin 9 

sin 9 cos 9 



e t/2 
t/2 



e -t/2 



1 « 
1 



n=0 



e*/ 2 (cos^ + nsin( 
e*/ 2 sintf 



n=0 



— log (e* (cos 9 + n sin 9) 2 + e* sin 2 6> ) 

— (t + log(l + n sin 20 + n 2 sin 2 0)) 
an 



n=0 



n=0 



sin 26> + In sin 2 # 



n=0 



1 + nsin(26>) + n 2 sin 2 9 
= sin(2#) 

Thus the critical points must satisfy sin (29) = 0, which means 9 = 
or 9 = 7r/2. Similarly we compute 

3 t/2 



(9a 



a=0 



d_ 

da 
d_ 

da 
d_ 



cos 9 — sin 9 



sin # 

e t/2 



COS 6* 

--t/2 



-t/2 



1 U 
1 



1 U 
1 



e'/ 2 (cosa + wsina) * 

g-t/2 g j n a s|e 



cos ct -sma 
sin a cos a 



a=0 



cos a 
sin a 



a=0 



2u 



log ^e'(cosa + wsina) 2 + e *sin 2 aj 



a=0 



and so the critical points of (p are (u,9) = (0,0) and (u,9) = (0, 7r/2) 
(in .fTAiV-coordinates); note however that these are simply the points 
based at ie l G H, pointing to oo and to — these points also have 
coordinates (0, t, 0) and (0, £, 7r/2) in iVAfr-coordinates. 

Now recall that 

K { ^(x,y,9) = I [ X {x-n)F L {2\9 -a\)e {is -^ n ^dadn 

We wish to apply the method of stationary phase to the integral over 
a,n. For this we will also need to compute the determinant \d 2 ip\ of 
the Hessian at the critical points (0,£, 0) and (0,£, 7r/2). 



n=0 



- sma 

cos a 
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First we consider the point (0, t, 0), based on the imaginary axis and 
pointing towards oo. Since <p is left-iV- invariant for vectors pointing 



towards oo, the second derivative 



dri- 



ft, and thus the second 



a=0 



derivative along a is irrelevant; indeed, the determinant \d <p\ of the 

8 2 <P 2 



dnda 



Hessian is simply 

This derivative is readily computed using ([3]); taking derivatives first 
with respect to n (at n — ft) and then with respect to a gives 



d 2 <p 



dadn 



(0,t,0) 



d_ 

da 



sin(2a) 



a =o e * cos2 a + e * s i n a 



2cos(2a) 



sin 2 (2a)(e-*-e^ 



e* cos 2 a + e~ l sin a 
2e"* 



(e* cos 2 a + e~* sin a) 



a=0 



and hence 
(4) 



(0,i,0) 



-e 
2 



At the point (0, t, 7r/2), on the other hand, we find that 
0fy 2cos(2a) sin 2 (2a) (e 



dadn 



t_ e t) 



(0,t,7r/2) 



e* cos 2 a + e * sin 2 a 



-2e* 



(e*cos 2 a + e *sin 2 o;) 2 



Here the derivatives along N are no longer identically zero, and so we 
must also calculate 



d 2 <p 



da 2 



(o,W 2 ) 



d 2 <p 



dn 2 



(0,*,7r/2) 



<9 2 



log(e' cos z a + e * sin 2 a 



da 2 

d ( -2 sinh(t) sin(2a) 
da \ e* cos 2 a + e~* sin 2 a 

4sinh(t)e' = 2e 2 ' - 2 
d 2 



log(e-*(n 2 + 1)) 



dn 2 
d ( 2n 



n=0 



dn \ n 2 + 1 



n=0 



a=iv/2 



a=n/2 
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and thus 

|<9Vr V2 = |4e 2 '-4-4e 2 *|- 1 / 2 

(0,t,7r/2) 

\ 

Thus, recalling that </?(0, t, 0) = £ and (p(0,t,ir/2) = —t, we find by 
the method of stationary phase (eg., |Hor03t Theorem 7.7.5]) 



I I x{x-n)F L (2\d-a\)e {ls ~^ n ^ a) dadn 

Jn Ja&T 



V\dWo,t,o)\ 

+ e^^xWFdn - 26) + L ,(s~ 2 ) 

v/|9V(o,t,f)| 

= 2 -^e^ t X {x)F L {26) + ^Lle^-^ X (x)F L (n - 29) + L)t (s 

= (e ist F L {29) + e-* st F L (n - 26)) + L , t ( S - 2 ) 

The error term depends on 2 derivatives of Fl — which are bounded by 
0(L 3 ) — and of x, which are fixed. The dependence on t is bounded 
by a fixed number of derivatives of ip(-, t, •), and powers of \dip(-, t, 
away from the critical points, along with the factor e*/ 2 ; each of which 
are controlled by powers of e'*L Thus since L < r ]°^ C an d e'*' < r)j C ', 

we can reduce the error term to 0(s~ 2 ) times a fixed power of ; a 
crude bookkeeping of the various remainders shows that the estimate 
O(s~ 2 r* 00 ^ C ) is more than sufficient. □ 

Now we continue by performing the integration over the spectrum: 
Proposition 2. We have 

K U (x,t,9) = nx(x)\F L (29) + F L (7r-26)}e t / 2 g(t) + O(r] 00/c - 1 ) 



Proof: Proposition [TJ shows that 
k u (x, t, 9) = n X (x)e t/2 J h{s) tanh(vrs) (e ist F L (29) + e- ist F L (n - 29)) ds+0{r l ™ /C ~ 
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since by (j2J) the integral of the error term O(r^ 00 ^ C s 2 ) certainly satis- 
fies 



100/C 



y s l h(s)t&r\h.(irs)ds < rj 00 ^ 1 



Moreover, since 1 > tanh(7rs) > 1 — 2e 2ws , and for |s| < ~Tj we have 
|/i(s)| < |rj — s|~ 3 by the rapid decay of h, we estimate 



g(t) - / e lst h(s) tanh(ns)ds 



< 



< 



< 



e ist h(s)[l - tanh(7rs)]ds 

— tanh(7rs)](is 

/ h(s)ds + / /i(s) • (1 — tanh(7i\< 
J\s\<U- J\s\>^4- 



|r> — s\ ds + 



sl>-^- 



h{s)e- z ™ds 



h(s) 



< 



M>4 



100/C-l^ 



which is readily absorbed into the error term 0(r 

Similarly for the e~ lst FL(^ — 28) term, noting that g is even so that 



-ist 



h(s)ds = g(-t) = g(t) 



and we are done. □ 



Corollary 2. Assume C is sufficiently large (eg. C > 200 is suffi- 
cient). Then for \t\ < 2C~ l logrj and |2#|, \n — 28 \ > r • 5 ^ C ; we have 



r ■ 



■3/C 



Proof: Plug into Proposition [2] the estimate 

sin 

-5/c r > „io/c 



whenever \28\ > and L > r ■ , and similarly for i^(7T — 20). 

Since x an d fl^) are uniformly bounded and e*/ 2 < r 1 - ' , the main 



terms of the k^j asymptotic from Proposition [2] are Oir- 4 ^ c ). As long 

as C is large, the error term 0(r* C 1 ) is also small enough and the 
Corollary holds. □ 
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4. Localization and the Proof of Theorem CD 

Recall we have a closed geodesic £ C M, and a point p G £. We 
lift M to the upper-half plane in such a way that p is mapped to the 
origin, and £ is mapped to the imaginary axis {iy : y > 0}. We will 
often abuse notation and also use £ to refer to the geodesic path in 
S*M. 

We consider the "collars" Uj C PSL(2R) ^ NAK defined by 

{(M,0) : \x\ < 1, \t\ < CT 1 log TV, |20| < r7 5/c } 
{(x,i,0) : |x| < 1, \t\ < CHogrj, |vr - 26\ < rJ 5/C } 

jjnp u ^down 

for C sufficiently large. 

We show that the projection of Uj to S*M is one-to-one, with the 
obvious exceptions of the periodicity of the geodesic £: 

Lemma 2. Let Uj as above. Then 

(u r u^)nVcA 

In other words, if g\,g% E Uj project to the same point g G S*M, 
then gi G T.g 2 and thus gig^ 1 G T; the lemma then guarantees that 
gi = a>ki(£)-g2, where Z(£) is the length of the geodesic £, and the element 
of T in question is simply wrapping around the geodesic £ exactly k 
times. 

Proof: We give the argument for Uj p , the argument for ^ down is 
completely analogous (as is the argument that the projections of U^ p 
and Uf own are disjoint). 

First, we note that Uf C B-(2)AB + (rJ 4/C ), where B~(2) is a 2- 
ball along the stable direction iV, while B + (r~ i ^ c ) is a small r~ 4 ^ c ball 
in the unstable direction iV, and A is the diagonal subgroup. Since in 
Uj the diagonal component is at most C~ l log rj and the fC-component 
is very small < r • , this means that up to minor adjustments to 
the constants, the i^-component can be replaced by NAN coordinates 
that are simpler to manipulate, since we intend to conjugate by large 
elements of A. 

Thus, suppose we have 

n{uJ)a t n{u + ) G Yn{vJ)a s n{y + ) 

where the stable components \u-\, < 2, the unstable components 
< r~ 4 / c , and |t|,|s| < 2C~ l \ogrj. We then conjugate by 



u7 = 

rrdown 
3 

Uj = 
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Ofc.j(5) G T, for a large integer /c of size ^fey- logr,-. This implies 

a t n(e- km+t u^)n(e km u + ) G raXe~ W(0+ V-)n(e w ®v + ) 
and moreover since ann G T we can write 

a T n(e- M(f)+t M_)n(e M(5) M + ) G ra CT n( e - fci(€)+s t;_)n(e M(f) t; + ) 
for |r|, | cr | < = 0(1). But note that 

| e~ kl ® +t u- 1 , | e- kl ^ +s v- 1 < rf C r 2 ' C < rJ 1/C 
\e kl ^u + \,\e kl ^v + \ < r] /C r- A/C < rJ 1/C 

are all small, and so this element of T is rj 1 -close (on the right) to 
the geodesic £. This implies that it is in fact on the geodesic itself, 
once Tj is large, since £ has a (right-)neighborhood in S*M that does 
not intersect any elements of T other than those in A (corresponding 
to the periodicity of £). □ 

Proposition 3. Let U D £ be an open set in S*M containing the 
geodesic £. Then for rj sufficiently large, 

\\ R u\\h(u) £ 1/logrj 
Proof: The neighborhood U contains a subset V C U of the form 

for some rj > 0, where we set 6* = to be a direction along the geo- 
desic contained in U. In this parameterization V is contained in the 
projection Uj p of Uj p to S*M, and disjoint from the projection JJ^ own 
of f/f™ 11 . Also note that since the contributions of K^j from points 
outside Uj p and f/j 10 ™ 11 are each bounded by rJ 3 ^ C by Corollary [21 and 

since there are 0(r • ) such points in the support of k^j — contained 
in the ball of radius C~ l log rj — mapping to each point of S*M (see 
eg. |LP82j ). these contributions will be negligible with respect to the 

final 1/logTj bound. Thus we must estimate the projection of k^j 
to V. 

Recall that the point 

t + za a a f et/2 ° "\ i f et/2 ° "\ a et A , t r u 



LOG-SCALE QUASIMODES THAT DO NOT EQUIDISTRIBUTE 19 

and so for a point given, in local ANK coordinates around p G £, by 

. fe T / 2 \ fl y\ (cos6 -sin0\ T _ 

the matrix _ /9 , • n n 6 k we nave 

V e / V V \ sm ^ cos ^ / 




\ (I y\ (cosd -sin0\\ 
e- T ' 2 ) \ I J \sin6 cos 6 J J 



trivially estimating the contributions from K^j outside Uj p by 0(r ■ ), 
as discussed above. Therefore we find by Proposition [2] and the defini- 
tion of g(t) in section [2721 



l%illW)-(l + 0(r7 1/C )) 

> r f [ ( J2 K U (ye kl ^ + \kl(0 + T,6)) drdydO 
J T =o J\v\<r,J\o\<r-*/a \\km<c~no grj J 

> fV^yr ( T x (ye kl ^ + ^ m0+ A dydr 



since J. & , <r -5/o F 2 {2d)d0 ~ 1 because L > r- , the symmetry in \ 

allows us to integrate only over the positive values of y at the expense 
of a constant factor of 2, and restricting the range of k decreases the 
overall value of the integrals, while allowing us to replace x(i^7~(^(0 + 
t)/ cosh( c ^^) +T ) ) > 1 in the definition of g{kl{^) + r)on this range. 

Note furthermore that since rj is a whole multiple of 2-7r/(£), we have 
cos (rj (£;/(£) + r)) = cos(r\,r), and can be pulled out of the sum over 
k; here we use our assumption that {r\,} is a sequence of frequencies 
tailored to the geodesic £. 

Moreover, further restricting the range of y and noting that x = 1 on 
[— 5) |]j an d that the r-dependence inside the sum can be abandoned 
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at the expense of another constant depending only on /(£), we find 

ll^llW)-(i + o(r7 1/c ; 



2 



•' ' \-^C-ilogri<fcI(e)<-Iog(2v) / 

(e- lo ^/ 2 ) 2 dy 



> 



1 


log 


; 2 fo) 




i 


log 


; 2 fo) 




i 


log 


; 2 fo) 




i 


log 


5 2 (^) 




1 



> X / ( -log(22/)/2^ 



-1/2C 

-vac 2y 



> r 2/ \ ( lo s ^ + lo s r i) £ 1 / c lo s r i 



l 

log^(r,) v "°'' ' 2C 
where we have used the fact that ye kl ^> < | when kl(£) < — log(2y), 

„ 1/2(7 

and have restricted the integral to the range y >r- to ensure that 
— log(2y) < iC _1 logrj in changing the sum over k. 

Since ||/c^j||L2(t/) > H%,ilU 2 (y)> the result follows. □ 
Proposition 4. VFe have for all rj sufficiently large 

\\ R u\\h(s*M) £ l/!og r j 
Combined with Proposition^ this implies 

|i- 1 12 1 1 - II 2 >1/1 "> 1 1 - II 2 

W^jW L 2 (S*M) — \ \ K Z,j\\L 2 (U) ~ l l l0 S r j <j \ \ K t,j\\L 2 (S*M) 

for any neighborhood U o/£; thus any quantum limit /i of the probability 
measures ^j(f) = p^-ypr j s , M f\R^j\ 2 dVol satisfies > 0. 

Proof: Again, since k^j < r~ 3 ^ c off of Uj by Corollary [21 and the 

projection of supp(^j) to S*M is 0(r - )-to-one, we may restrict our 
attention to the projection of Uj to S"*M; i.e. 

— 1/C 

where t/J 13 and U^ own are the respective projections to 5"*M of Uj p and 
jjdawn^ an( j we can evaluate the two norms on the right separately; since 
they are essentially identical, we will do the Uj p term. Note that 

1 y\ fcosO -sin0\ , , i/c la] . -5/c 
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We follow the calculations in the proof of Proposition |3j this time, 
though, we trivially estimate all % terms, as well as the cos / cosh term, 
and L Fl(29)d9, by 1. We recall however that x is supported on [—1, 1], 
and so 



-1(0 

< 



< 



r=0 J\y\<r 1/C J\6\<r 



-5/C 



£ K U {ye km+ \ kl(0 +T,e)\ drdyde 



\W(0<C- 1 logr j 



log ( r j) Jt=0 J y =o 



kl(0+r\Mkl(0+r) 



dydr 



< 



.i/o 



log (rjj jy= r . 

1 



1/C 



< 



< 



+ 



1 



Jk\l(0<C-Ho srj 

-C" 1 logrj<kl(£)<-\ogy 
-I/O 



log (rj) J y=0 

1/C 



log (r 



3 ) \Jy=r 



-1/C 



e^ km 

\k\l(0<C~i lo grj 
r~ 1/C 

e- logy dy+ I' r) /C dy 

Jy=0 



dy 



— —(C-Mogr,- + 1) < 1/logr,- 

log [Tj) 



as required. The [/? own term is identical, using the analogous trivial 
estimate of f $ Fl(ir - 26)d6 < 1. □ 

Finally, we must verify that k^j is indeed an O (^^r j-quasimode. 

From the above calculations in Proposition HI it is most convenient to 
do this via the estimates for k^j. 

Proposition 5. We have 



< -^-llfe 

2 ~ l0 S r i 



and thus {k^j} is a sequence of O ( j^-p- j-quasimodes. 
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Proof: Since | |l 2 (m) = I \l 2 (s*m), we can replace the right- 



hand side with 
pute on H 



T^F-\\ K aj\\L2{s*M) 



1 

4 



A + t + ] /s w 



. For the left-hand side, we first com- 



s 2 )kg jSh(s) tanh(7rs)o?s 



Since the operators R n and R~ n defining the microlocal lift in Lemma [T] 
act in each irreducible representation of spectral parameter s, and (A+ 
| + rj) also commutes with T-translations on the left, we can replace 
the left-hand side of Proposition [5] with 



L 2 (M) 



E 
E 

7er 



s 2 )k^ ji'j-)sh(s) tanli7rs(is 



s 2 )k^](7-)s/i(s) tanhvrsc/s 



L 2 (M) 



L 2 (S*M) 



and then we may perform all calculations in terms of the microlocal 

(s) 

lifts k£ L where we have already estimated the stationary phase asymp- 
totics. 

Indeed, by Proposition [TJ and following the proof of Proposition |2l 



r| — s 2 )/^ s j(x, i, 6)sh(s) ta,nh(7rs)ds 



t/2 



r 2 - s 2 )h(s) e tst F L (26) + e~ tst F L (-K - 26) )ds + 0(r 



nx{x)e 



n X (x)e^(F L (2B) + F L (n - 26)) ■ (r 2 g(t) + g"(t)) + 0( 



,100/C-1> 



[r 



100/C-lN 



and it remains to estimate the function r 2 g(t) + g"(t). 
Recall that 



C 



cos(trj 



X 



c 



log r j 



where K« 



21o. 



41 °^ cosh(^| 
= cos(tr i ) • KH(Kt) 

c and := |x(20/ c osh<. Thus 



-Tj sin(tr j ) • K H'(Kt) + 0(K 6 H"(Kt)) 

1 



—(rjK) ■ sin{trj)KH'{Kt) + O 



log 3 r j 
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Now sin(trj) and H' are uniformly bounded in rj, so that the calcu- 
lations of Proposition 0] — where g(t) was estimated trivially in the sup- 
port of H — may be applied equally well to the function sin (rjt)KH'(Kt) 
in place of g(t) = cos(rjt)K H (Kt) , and thus we arrive at 



A 



1 2 

7 + r i 
4 3 



k. 



< 



< 



< 



,j. 2 



log Tj 
j. 2 

log 2 r j 



log Tj 

\tc II 2 

\ K i,j\\L 2 (S*M) 

\b II 2 
l^UW L 2 (M) 



as required, using the estimate \\k^ 



II 2 

,j\\2 



1/logrj from Proposition H] 



and the fact that H/^jHi^r^M) = | \l 2 (m)- D 



Combining Propositions El HI and |5l we have now proved Theorem [T] 



for 



a 

log r 



-quasimodes, with a sufficiently large constant C = 0(C); 

in the next section we will describe how to build on this case to reduce 
the constant from a large C to arbitrarily small e > 0. 



5. Finer-scale Quasimodes 

We now discuss how to use these {j^-j- quasimodes k^j, that con- 
centrate at least Si > of their mass on the geodesic £, to construct 
^i^j -quasimodes that concentrate positive mass on £, for any e > 0. 

Naturally, in agreement with Conjecture El our lower bound on the 
mass concentration will shrink as e — > 0. 

The idea is to start with an arbitrary o ^j^r j-quasimode ipj near 
rj — for example, an eigenfunction with spectral parameter near r 3 — 
whose microlocal lift tyj we may assume does not concentrate mass 
on £, and perturb it by fak^j. The sum ipj + #2%j will be an j^r- 
quasimode if 5 2 is small enough, and \&j + ^2%j will still concentrate 
5 2 Si of its mass on £. 

In order for this to work, ipj must have approximate eigenvalue near 
Tj] precisely, within K/logrj, for some fixed K independent of e. By 
the estimates of Berard |Ber77| for the remainder term in Weyl's Law, 
we know that there exists a K depending only on the surface M, such 
that an interval of width K/logr around r must contain at least one 
eigenvalue. So we may indeed take ipj to be such an eigenfunction, 



24 SHIMON BROOKS 

with spectral parameter fj satisfying 

< fj — Tj < Kj log Tj 

and set 82 = e/(2K + C). Naturally, we normalize H^jlh — ll^jlb- 

If (hypothetically) the microlocal lifts ^fj of these ipj were to concen- 
trate > ^8281 of its mass on £, then there is nothing to prove; these ipj 
are eigenfunctions, and thus in particular are j^r- quasimodes for any 

e > 0. Otherwise, we define kj = ipj + 5 2 &f j, whereby its microlocal 
lift is 

>\, = *i + 'V.'i.., 

Since by the triangle inequality, for any sufficiently small neighbor- 
hood U D £ in S*M and j sufficiently large, we have 

||^il|L 2 (C7) > \\82^,j\\L 2 {U) - \ \^j\\L 2 (U) > -^8 2 8i 

we see that kj concentrates a positive proportion of its mass on the 
geodesic £. It remains to show that kj is indeed an quasimode. 

But this is an immediate consequence of the fact that we diluted k^j 
by a sufficiently large factor. Observe that since ipj is an eigenfunction 
of spectral parameter fj we have 

||(A+(i+fJ))^|| 2 < ||(A+(i+f|M-|| 2 + 5 2 ||(A + (i + f|))^| 2 

1 



< 52\\(A+q+f]))k u \\ 2 

^(||(A + (i+r|)fe|i 2 + ir|-f|i-|fe|i 2 



< 



< 82 + 



log Tj log Tj 

e 



J 1 

log r j 

since 8 2 = e/(2if + C"), and < fj — r 3 - < K/logr.,, completing the 
proof of Theorem [TJ □ 
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